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1 Introduction 



The ratio of Cabibbo-Kobayshi-Maskawa matrix [T] elements |Kd|/|Ks| has been measured 
in B mixing processes, where it is probed in box diagrams through the ratio of and 
Bg mixing frequencies [2]-|5]. The ratio of these matrix elements has also been measured 
using the ratio of branching fractions of 6 — )■ s'j and b ^ d'j decays, where radiative 
penguin diagrams mediate the transition |6]-{8]. These measurements of | VJ;d|/| Ks| are 
consistent, within the (dominant) ~10% uncertainty on the determination from radiative 
decays. The decays b — )• s^'^fi~ and b — )■ dfi'^fi~ offer an alternative way of measuring 
|Kd|/|K;s| which is sensitive to different classes of operators than the radiative decay 
modes |9]. These 6—)- {s,d)^~^^~ transitions are flavour-changing neutral current processes 
which are forbidden at tree level in the Standard Model (SM). In the SM, the branching 
fractions for 6— >■ d£~^i~ transitions are suppressed relative to 6— )■ si~^i~ processes by the 
ratio |VtdP/|KsP- This suppression does not necessarily apply to models beyond the SM, 
and B'^ — )■ Tc^fi~^ij,~ decay^ may be more sensitive to the effect of new particles than 
B^ — )■ K^n^n^ decays. In the SM, the ratio of branching fractions for these exclusive 
modes 

R = B{B+ 7r+/i+;U-) / B{B+ K^fi+fi-) (1) 

is given by i? = F^/^, where V = | Kdl/I Ks] and / is the ratio of the relevant form factors 
and Wilson coefficients, integrated over the relevant phase space. A difference between 
the measured value of R and ^^/^ would indicate a deviation from the minimal flavour 



models 12 



violation hypothesis 10 ,11 , and would rule out certain approximate flavour symmetry 



No 6—)- di~^£~ transitions have previously been detected, and the observation of the 
B^ —7- Ti^fi^fi" decay would therefore be the first time such a process has been seen. The 
predicted SM branching fraction for Tx^ji^fi" is (2.0 ± 0.2)xl0~^ [13]. The most 



stringent limit to date is B{B~^ n'^^^fi ) < 6.9 x 10 ^ at 90% confidence level 14 . The 



analogous 6— >■ si~^i decay, B^ ^ K^n^jj, , has been observed with a branching fraction 



of (4.36 ± 0.15 ± 0.18) X 10"^ 15 



This paper describes the search for the B^ — > decay using pp collision 

data, corresponding to an integrated luminosity of 1.0 fb~^, collected with the LHCb 
detector. The B^ — )■ 7r"'"/x^/i~ branching fraction is measured with respect to that of 
B^ ^ J/ip{-^ /i+/i~)ii'+, and the ratio of B^ ^ 7r"'"/x"'"/x~ and B^ ^ ix'^ ii~ branching 
fractions is also determined. 



The LHCb detector 16 is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < 77 < 5. The experiment is designed for the study of particles containing b 
or c quarks. The apparatus includes a high precision tracking system, consisting of a silicon- 
strip vertex detector surrounding the pp interaction region, and a large-area silicon-strip 
detector located upstream of a dipole magnet. The dipole magnet has a bending power 
of about 4 Tm. Three stations of silicon-strip detectors and straw drift-tubes are placed 
downstream of the magnet. The combined tracking system has a momentum resolution 

^Charge conjugation is implicit throughout this paper. 
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Ap/p that varies from 0.4% at momenta of 5GeV/c, to 0.6% at lOOGeV/c. The tracking 
system gives an impact parameter resolution of 20 |J.m for tracks with a high transverse 
momentum (px)- Charged hadrons are identified using two ring- imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a system composed of alternating layers 
of iron and either multi-wire proportional chambers or triple gaseous electron multipliers. 

In the present analysis, events are first required to have passed a hardware trigger 
which selects high-px single muons or dimuons. In the first stage of the subsequent software 
trigger, a single high impact parameter and high-px track is required. In the second stage 
of the software trigger, events are reconstructed and then selected for storage based on 



either the (partially) reconstructed B candidate or the dimuon candidate 17,18 



To produce simulated samples of signal and background decays, pp collisions are 



generated using Pythia 6.4 19 with a specific LHCb configuration 20 . Decays of 



hadronic particles are described by the EvtGen package [2T| in which final state radiation 



is generated using Photos 22 . The interaction of the generated particles with the 
detector and the detector response are implemented using the Geant4 toolkit 23 , as 



described in Ref. 24 



The small branching fractions of the — > Tr+yU^/i" and B^ — )■ K~^fi'^fi~ signal decays 
necessitate good control of the backgrounds and the use of suitably constrained models to 
fit the invariant-mass distributions. The decay B^ ^ J/ip[-^ ii~)K^ (hereafter denoted 
B^ — 7- J/ipK^) is used to extract both the shape of the signal mass peaks and, in the 
B^ —7- n^fi^fi'^ case, the invariant mass distribution of the misidentified B^ — ?► K^fi^fi~ 
events. These misidentified B^ — t- ^~ events form the main residual background 

after the application of the selection requirements. 



2 Event selection 

The B^ —J- -n^ jj,^ and B^ — )• i^^/i+zi^ candidates are selected by combining pairs of 
oppositely charged muons with a charged pion or kaon. The selection includes requirements 
on the impact parameters of the final-state particles and B candidate, the vertex quality 
and displacement of the B candidate, particle identification (PID) requirements on the 
muons and a requirement that the B candidate momentum vector points to one of the 
primary vertices in the event. The rate of events containing more than one reconstructed 
candidate is 1 in ~20,000 for B^ — )■ J/ipK^ . No restriction is therefore placed on the 
number of candidates per event. 

The pion identification requirements select a sample of pions with an efficiency of ~70% 
and a kaon rejection of 99%. The kaon identification requirements allow the selection of a 
mutually exclusive sample with similar efficiencies. The muon identification requirements 
have an efficiency of ~80%, with a pion rejection of ~99.5%. The PID requirements have 
a momentum dependent efficiency which is measured from data, in bins of momentum, 
pseudorapidity and track multiplicity. The efficiency of the hadron PID requirements is 
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measured from a sample of D*^ — )■ [D^ — i^'~7r^)7r^ candidates that allows the hadrons 
to be unambiguously identified based on their kinematics. The muon PID efficiencies are 
measured using t- J/ipK^ candidates, using a tag and probe method. 

The J/ip and i^i^S) resonances, where J/ip,ip(2S) — )■ are excluded using a veto 

on the dimuon mass. This veto has a total width of 200 (150) MeV/c^ around the nominal 



J/i/j {ip{2S)) mass 25 , and takes into account the radiative tail of these decays. Several 
other backgrounds are considered: combinatorial backgrounds, where the particles selected 
do not originate from a single decay; peaking backgrounds, where a single decay is selected 
but with one or more particles misidentified; and partially reconstructed backgrounds, 
where one or more final-state particles from a B decay are not reconstructed. These 
backgrounds are each described below. 



2.1 Combinatorial backgrounds 



A boosted decision tree (BDT) |26| which employs the AdaBoost algorithm 27 is used to 
separate signal candidates from the combinatorial background. Kinematic and geometric 
properties of the candidate and final state particles, candidate vertex quality and 
final state particle track quality are input variables to the BDT. 

The BDT is trained on a simulated B~^ — > ti^ signal sample, and a background 
sample taken from sidebands in the B^ — )■ 7r+/i+/i~ and B^ — )■ K^fi'^fi~ invariant 
mass distributions. These invariant masses are denoted M^+^+^- and M^+^+fM-, re- 
spectively. The background sample consists of 20% of the candidates with Mjr+^+^- or 
Mx+^i+fj.- > 5500 MeV/c^. This sample is not used for any of the subsequent analysis. 
Signal candidates are required to have a BDT output which exceeds a set value. This 
value is determined by simulating an ensemble of datasets with the expected signal and 
background yields, and choosing the cut value which gives the best statistical significance 
for the B'^ — > 7r+/i+/i~ signal yield. The same method is used to select the optimal set 
of PID requirements. The BDT output distribution for simulated B^ — t- 7r+|U+/i~ events 
and for mass sideband candidates is shown in Fig. [T} A cut on the BDT output > 0.325 
reduces the expected combinatorial background from 652 ± 11 to 9 ± 2 candidates in a 
±60 MeV/c^ window around the nominal B mass, while retaining 68% of signal events. 
Assuming the SM branching fraction and the single event sensitivity defined in Sect. |4| 
21 ± 3 i?+-> /T+yU+/i~ signal events are expected in the data sample. 



2.2 Peaking and partially reconstructed backgrounds 

Backgrounds from fully reconstructed B~^ decays with one or more misidentified parti- 
cles have a peaking mass structure. After applying the PID requirements, the fraction 
of B^^K^^^fi" candidates misidentified as B~^ ^ tt^ fi~ is 0.9%, giving a resid- 
ual background expectation of 6.2 ± 0.3 candidates. This expectation is computed by 
weighting B^ — )■ K~^fi'^fi~ candidates, isolated using a kaon PID requirement, accord- 
ing to the PID efficiency obtained from the D*~^ calibration sample. The only other 
decay found to give a significant peaking background in the search for B^ — > 7i^fi~^n~ is 
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BDT output 



Figure 1: BDT output distribution for simulated — )• 7r+/x+/i^ events (black solid 
line) and candidates taken from the mass sidebands in the data (red dotted line). Both 
distributions are normalised to unit area. The vertical line indicates the chosen cut value 
of 0.325. 



—7- TT+TT+TT^, where both a tt^ and a vr^ are misidentified as muons. For B^ — )■ K^fi^jj' 
decays, the only significant peaking background is B^ — t- K^n^n^ , which includes 
the contribution from D°(— K~^7[~)tt~^. The expected background levels from 

B+ TT+TT+TT- {B+^ K+TT+TT-) dccays are computed to be 0.39 ± 0.04 (1.56 ± 0.16) 
residual background candidates, using simulated events. 

Backgrounds from decays that have one or more final state particles which are not 
reconstructed have a mass below the nominal B mass, and do not extend into the signal 
window. However, in the B^ — > 7i^fi^fi~ case, these backgrounds overlap with the 
misidentified B^ — )■ K^fi'^fi" component described above, and must therefore be included 
in the fit. In the B^ — )■ K^fi^fi" case such partially reconstructed backgrounds are 
negligible. 

2.3 Control channels 

The B^ —7- J/ip and B^ — decay candidates are isolated by replacing the pion 
FID criteria with a requirement to select kaons. In addition, instead of the dimuon mass 
vetoes described above, the B^ — )■ J/ip candidates are required to have dimuon mass 
within ±50 MeV/c^ of the nominal J/ip mass (the J/ip mass resolution is 14.5 MeV/c^). 
The remainder of the selection is the same as that used for B^ — )■ 7r+/i+/i~ . This minimises 
the systematic uncertainty on the ratio of branching fractions, although the selection 
is considerably tighter than that which would give the lowest statistical uncertainty 
on the B'^ — )■ K^fi'^fi^ event yield. The B^ — )• {J/ip — )■ ji^ ji^)-n^ candidates (denoted 
B^ —7- J/ifj-K^)^ which are discussed below, are selected using the same BDT, the pion FID 
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criteria, and the above window on the dimuon invariant mass. There is no significant 
peaking background for — )■ J/tp decays. For — >■ J/i(j 7r~^ decays the only significant 
peaking background is misidentified — )■ J/ip events. 

3 Signal yield determination 

The B^ —7- 7r+yU^/i~ , B^ — )■ K~^fi'^fi~ and B'^ — )■ J/tp yields are determined from a 
simultaneous unbinned maximum likelihood fit to four invariant mass distributions which 
contain: 

1. Reconstructed S"*"— )■ J/ipK~^ candidates; 

2. Reconstructed )► J/ipK^ candidates, with the kaon attributed to have the pion 
mass; 

3. Reconstructed 5+— > -n^ji^ji^ candidates; and 

4. Reconstructed B^ — )■ /^+/i+yU~ candidates. 

The signal probability density functions (PDFs) for the B^ — > vr^/i^/i", B^ — )■ K~^fi'^^~, 
and B^ — )■ J/ipK^ decay modes are modelled with the sum of two Gaussian functions. 
The PDFs for all of these decay modes share the same mean, widths and fraction of 
the total PDF between the two Gaussians. The B^ TT~^fi~^fi~ PDF is adjusted for the 
difference between the widths of the B^ — > n^fi'^fi^ and B^ — t- J/ipK^ distributions, which 
is observed to be at the percent level in simulation. The peaking backgrounds described 
in Sect. 2.2 are taken into account in the fit by including PDFs with shapes determined 
from simulation. The combinatorial backgrounds are modelled with a single exponential 
PDF, with the exponent allowed to vary independently for each distribution. The partially 
reconstructed candidates are modelled using a PDF consisting of an exponential distribution 
cut-off at a threshold mass, with the transition smeared by the experimental resolution. 
The shape parameters are again allowed to vary independently for each distribution. The 
misidentified i?"*"— )■ J/ipK^ candidates are modelled with a Crystal Ball function |28], as 
it describes the shape well. In order to describe the relevant background components for 
B^ —7- 'K^fi~^fi~, the fit is performed in the mass range 4900 < M^+^+^- < 7000 MeV/ c^. To 
avoid fitting the partially reconstructed background for B^ — )■ K^fi^fi" , which is irrelevant 
for the analysis, the fit is performed in the mass range 5170 < Mk+^+^- < 7000 MeV/c^. 

3.1 Reconstructed J/il^K^ candidates 

The reconstructed B^ — t- J/ipK^ candidates are shown in the M^+^+^- distribution 
in Fig. |2fa). The fitted 5+^ J/^K^ yield is 106,230 ± 330. This large event yield 
determines the lineshape for the B^ — )■ 7c'^fi~^fi~ and B~^ — )■ K~^fi~^fi~ signal distributions, 
and provides the normalisation for the B~^ 7r~^ fi~^fi~ branching fraction. 
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Figure 2: Invariant mass distribution for 5+—)- J/ipK^ candidates under the (a) K^fi~^fi~ 
and (b) n'^fi^fi" mass hypotheses with the fit projections overlaid. In the legend, "part, 
reco" refers to partially reconstructed background. The fit models are described in the 
text. 



3.2 Reconstructed B'^ J/iIj K'^ candidates with the pion mass 
hypothesis 

The — )■ J/ipK^ candidates reconstructed under the pion mass hypothesis provide 
the lineshape for the misidentified — )■ K~^fi'^fi~ candidates that are a background 
to the B^ —7- TT+yU+yU^ signal. The equivalent background from B^ — )■ 7r+/x+/x^ in the 
B^ —7- K^fi^fi~ sample is negligible. 

The PID requirements used in the selection have a momentum dependent efficiency 
and therefore change the mass distribution of any backgrounds with candidates that have 
misidentified particles. In order to correct for this effect, the B~^ ^ J/ipK^ candidates are 
reweighted according to the PID efficiencies derived from data, as described in Sect. 2.2. 
This adjusts the 5+— t- J/ipK^ invariant mass distribution to remove the effect of the kaon 
PID requirement used to isolate B^ ^ J/ipK^, and to reproduce the effect of the pion 
PID requirement used to isolate B^ ^ 7r+/i+/i~. In addition, there is a difference in the 
lineshapes of the B^ — )• J /if) and B^ — )■ K^fi^fi^ invariant mass distributions under the 
pion mass hypothesis. This effect arises from the differences between the two decay modes' 
dimuon energy and hadron momentum spectra, and is therefore corrected by reweighting 
S^— )■ J/ipK^ candidates in terms of these variables. The M^+^+^- distribution after both 
weighting procedures have been applied is shown in Fig. [2]^b). 

3.3 Reconstructed — )- Tr+^t+^i," and B^ — )- ^~ candi- 
dates 

The yield of misidentified B^ K^^'^^~ candidates in the B^ vr+yU+yU^ invariant mass 
distribution is constrained to the expectation given in Sect. 2.2. Performing the fit without 
this constraint gives a yield of 5.6 ± 6.4 misidentified i?"*"— )• K~^fj,'^fj,~ candidates. The 
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Figure 3: Invariant mass distribution of 7r+/i+/i^ candidates with the fit projection 

overlaid (a) in the fuU mass range and (b) in the region around the B mass. In the legend, 
"part, reco." and "combinatorial" refer to partially reconstructed and combinatorial 
backgrounds respectively. The discontinuity at 5500 MeV/ is due to the removal of data 
used for training the BDT. 




Figure 4: Invariant mass distribution of B~^ — )■ K^fi'^fi~ candidates with the fit projection 
overlaid (a) in the full mass range and (b) in the region around the B mass. In the legend, 
"combinatorial" refers to the combinatorial background. 



yields for the peaking background components are constrained to the expectations given in 
Sect. 2.2. For both the M^+^+^- and Mk+^+^- distributions, the exponential PDF used 
to model the combinatorial background has a step in the normalisation at 5500 MeV/c^ 
to account for the data used for training the BDT. 

The M^+^+^- and M^+^+^- distributions are shown in Figs [3] and |4| respectively. The 
fit gives a 5+ — > 7r+/i+yU^ signal yield of 25.3 Ig'J, and a 5+ — t- K^fi^fi^ signal yield of 
553 _25. 
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Figure 5: Invariant mass distribution of — )■ J/ip-K^ candidates with the fit projection 
overlaid. In the legend, "part, reco." and "combinatorial" refer to partially reconstructed 
and combinatorial backgrounds respectively. The fit model is described in the text. 



3.4 Cross check of the fit procedure 

The fit procedure was cross-checked on — t- J/ipir^ decays, accounting for the background 
from B^ — )■ J /if) decays. The resulting fit is shown in Fig. [5j The shape of the combined 
B^ — )■ J/ipTc^ and B^ — )• J/ipK^ mass distribution is well reproduced. The B~^ — )■ J/ipK~^ 
yield is not constrained in this fit. The fitted yield of 1024 ± 61 candidates is consistent 
with the expectation of 958 ± 31 (stat.) candidates. This expectation is again computed by 
weighting the -B"*"— )■ J/ipK^ candidates, which are isolated using a kaon FID requirement, 
according to the FID efficiency derived from D*^ — )■ {D^ — )■ K~tt^)tt^ events. 



4 Determination of branching fractions 

The B^ — 7r+/x+/x^ branching fraction is given by 

B{B TT fl fl ) = B+^7T+ ,1+ ti- (,2) 

= a ■ NB+^n+tJ.+ii- , (3) 

where B{X), Nx and ex are the branching fraction, the number of events and the 
total efficiency, respectively, for decay mode X, and a is the single event sensitiv- 
ity. The total efficiency includes reconstruction, trigger and selection efficiencies. The 
ratio eB+^J/^K+/^B+^■K+^l+^i- is determined to be 1.60 ± 0.01 using simulated events, 
where the uncertainty is due to the limited sizes of the simulated samples only. Other 
sources of systematic uncertainty are discussed in Sect. [5j The difference in efficien- 
cies between _B+ — )■ 7r+/i+/i~ and B^ — )■ J/ipK^ events is largely due to the mass vetoes 
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used to remove the charmonium resonances, and the different PID requirements. The 
J/i){-^ branching fraction is (6.02 ± 0.20) x 10"^ |25]. Together with the 

other quantities in Eq. 2, this gives a single event sensitivity of a = (9.1 ± 0.1)xlO^^°, 
where the uncertainty is due to the hmited sizes of the simulated samples only. 
The ratio of — t- TT^fi~^fi~ and K^fi'^fi~ branching fractions is given by 

Ns+^n+ti+iJ.- ^B+~^K+fi+ti- 
^ ' ^> 

where simulated events give eB+^x+^i+fi- /^B+^n+fi+fi- = 1-15 ± 0.01. 



5 Systematic uncertainties 

Two sources of systematic uncertainties are considered: those affecting the determination 
of the B~^ — > 7r+/i+/i~ and B~^ — )■ K~^fi'^^~ signal yields, and those affecting only the 
normalisation. 

Uncertainties in the shape parameters for the misidentified B^ — > K^fi^fi^ PDF in 
the fit are taken into account by including Gaussian constraints on their values. The 
most significant sources of uncertainty in the determination of these shape parameters 
arise from the procedure for correcting the B~^ — )■ J/ipK^ mass shape to match that of 
the B^ — )■ K'^n'^fi' decay, and the correction for the hadron PID requirements. The 
uncertainty on the B^ — TT~^fi~^fi~ yield determined with the fit takes these shape parameter 
uncertainties into account, and they are therefore included in the statistical rather than 
the systematic uncertainty. These uncertainties affect the B^ — t- 7r+|U+/i~ yield at below 
the one percent level. None of these effects give rise to any significant uncertainty for the 
B~^ — )■ K~^fi~^fi~ decay. 

Uncertainties on the two efficiency ratios ^B+^j/ipK+ / ^B+^■K+^l+^l- and 
^B+^K+^l+^l- |^B+^■K+^l+^l- affect the conversion of the B^ — )• 7r+yU+yU^ yield into a 
branching fraction, and the measurement of the ratio of branching fractions R. The 
largest systematic uncertainty on these efficiency ratios is the choice of form factors used 
to generate the simulated events. Using an alternative set of form factors changes the 
i?"*"— )■ Tr^/i+zi" efficiency by 3%, and this difference is taken as a systematic uncertainty. 
For the ratio of 5+— )■ TT^fi^fi~ and 5+— ?■ K^fi^fi~, the alternative form factors are used 
for both B^ 7r+|U+yU^ and B^ — )■ K^^^^^ , giving a systematic uncertainty of 1.7%. 
To estimate the uncertainty arising from the PID efficiency, the ratio of corrected yields 
between the -B"*"— ?• J/ipK^ and -B^— )■ J/ipir^ decay modes is measured, varying the PID 
requirements. The largest resulting difference with respect to the nominal value is 1.1%, 
which is taken as the systematic uncertainty. 

The systematic uncertainty arising from the knowledge of the trigger efficiency is 
determined using B^ — t- J/ip candidates in the data. Taking the events which pass 
the trigger independently of the B^ — )■ J/ip candidate, the fraction of these events 
which also pass the trigger based on the B^ — )■ J/il) candidate provides a determination 
of the trigger efficiency. The efficiency determined in this way is compared to that 
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Table 1: Summary of systematic uncertainties. 



Source 




^3iB+^K+^J.+fl-) y^'^) 



Form factors 

Trigger efficiency 

PID performance 

Data simulation differences 

Simulation sample size 

Total 



1.7 
1.4 
1.1 
0.4 
0.7 



calculated in simulated events using the same method, and the difference is taken as 
the systematic uncertainty. This gives a 1.4% uncertainty on eB+^j/ii,K+ /^B+^-K+iM+fM- and 

^B+^K+fi+^j.-^ I ^B+'^TT+fi+fj.- ■ 

For all decays under consideration, there are small differences between the distributions 
of some reconstructed quantities in the data and in the simulated events. These differences 
are assessed by comparing the distributions of data and simulated events for — )■ J/tp 
candidates. The simulation is corrected to match the data where it disagrees, and the 
resulting 0.4% difference between the raw and corrected ratio of — )■ J/ipK^ and 
B^ — )■ Tr+yU^/i" efficiencies is taken as a systematic uncertainty. The statistical uncertainty 
from the limited simulation sample size is 0.7%. When normalising to S"*"— )■ J/ipK^, the 
measured 5+— t- J/ipK^ and J/ip ^ ^^fi~ branching fractions contribute an uncertainty 
of 3.5% to the B^ — )■ n^fi^fi'^ branching fraction. The systematic uncertainties are 
summarised in Table [H 



6 Results and conclusion 

The statistical significance of the B^ — t- n^fi^fi^ signal observed in Fig. |3] is computed 
from the difference in the minimum log-likelihood between the signal-plus-background 
and background-only hypotheses. Both the statistical and systematic uncertainties on 
the shape parameters (which affect the significance) are taken into account. The fitted 
yield corresponds to an observation of the i?"*"— )■ 7r"'"/i"'"/i~ decay with 5.2 a significance. 
This is the first observation of a 6—)- di^i~ transition. Normalising the observed signal to 
the B^ ^ J/ipK~^ decay, using the single event sensitivity given in Sect. |4| the branching 
fraction of the B^ n^fi^fi^ decay is measured to be 

B{B+^n+li+fi-) = (2.3 ± 0.6 (stat.) ± 0.1 (syst.)) x 10"^ 



This is compatible with the SM expectation of (2.0 ± 0.2) x 10'*^ 13 



Taking the measured 5+— i^+yU+yU yield and tB+^K+ti+ii- /^B+^n+ti+fi-, the ratio of 
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and 5+ i^+zi+zi" 
g(g+^7r+/iV-) _ 



branching fractions is measured to be 
0.053 ± 0.014 (stat.) ± 0.001 (syst.) 



In order to extract |Vtd|/|Vts| from this ratio of branching fractions, the SM expectation 
for the ratio of — )• 7r"'"/x+/i~ and B^ — )■ K~^fi~^fi~ branching fractions is calculated using 
the EvtGen package 21 , which implements the calculation in Ref. 29 . This calculation 



has been updated with the expressions for Wilson coefficients and power corrections from 
Ref. 30 , and formulae for the dependence of these coefficients from Refs. 31 , 32 . 
Using this calculation, and form factors taken from Ref. [33] ("set 11"), the integrated 
ratio of form factors and Wilson coefficients is determined to be / = 0.87. Neglecting 
theoretical uncertainties, the measured ratio of B^ — )> n^fi'^fj," and B^ — )• K^fi^fi^ 
branching fractions then gives 



Kdl/IK, 



ts| 



0.266 ± 0.035 (stat.) ± 0.007 (syst.) 



which is compatible with previous determinations |5j-[8] . An additional uncertainty will arise 
from the knowledge of the form factors. As an estimate of the scale of this uncertainty, 
the "set IV" parameters available in Ref. 33 change the value of |Vtd|/|Ks| by 5.1%. 



This estimate is unlikely to cover a one sigma range on the form factor uncertainty, and 
does not take into account additional sources of uncertainty beyond the form factors. 
A full theoretical calculation taking into account such additional uncertainties, which 
also accurately determines the uncertainty on the ratio of form factors, would allow a 
determination of iKdl/lKsl with comparable precision to that from radiative penguin 
decays. 
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